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Detection of compact objects by means of gravitational lensing 

in binary systems 
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Abstract. We consider the gravitational magnification of light for binary systems containing two compact objects: 
white dwarfs, a white dwarf and a neutron star or a white dwarf and a black hole. Light curves of the flares of the 
white dwarf caused by this efi'ect were built in analytical approximations and by means of numerical calculations. 
We estimate the probability of the detection of these events in our Galaxy for different types of binaries and show 
that gravitational lensing provides a tool for detecting such systems. We propose to use the facilities of the Sloan 
Digital Sky Survey (SDSS) to search for these flares. It is possible to detect several dozen compact object pairs 
in such a programme over 5 years. This programme is apparently the best way to detect stellar mass black holes 
with open event horizons. 
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1. Introduction 

One of the most important manifestations of gravitational 
lensing is the visible variability of the astrophysical objects 
whose emission is affected by the gravitation field of the 
lens. This effect becomes observable over a reasonable time 
period, when the velocities of relative motions of the ob- 
server, the lens and the source are great enough, i.e. their 
mutual location changes rapidly. Many cases of this kind 
of variability have been examined: from light variation of 
distant quasars to outburst of stars caused by the influence 
of planets (see, for instance, Zakharov 1997). The bright- 
ness variation of the binary system components caused 
by gravitational lensing was been first considered by Ingel 
(1972, 1974) and Maeder (1973). They noted that repe- 
tition and small characteristic times of the effect make it 
one of the most accessible to detection and detailed study. 
It has also been shown (Maeder 1973) that the large am- 
plitudes of brightness variations (CK"1 — 0™5) can be ex- 
pected in binary systems consisting of compact objects 
- white dwarfs, neutron stars, black holes. In fact, only 
in these cases does the Einstein-Hvol'son radius turn out 
to be smaller than the lens star radius, and the image of 
the radiation source is practically not overlapped by it. 
The latter circumstance makes the gravitational lensing 
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influence in binary systems an effective means for the de- 
tection and detailed investigation of compact objects with 
the aid of photometric methods alone. Their permeabil- 
ity is by 2™ — 3™ higher (as compared to spectroscopy), 
and, therefore, the number of stars amenable to study is 
an order of magnitude larger. Here the light curve anal- 
ysis determines all the parameters of the binary system 
in full analogy to this task for eclipsing binary systems 
(with allowance made for its singularity) (see, for exam- 
ple, Goncharsky et al. 1985). 

A complete set of data for compact object binaries 
and their statistical analysis gives a possibility to test bi- 
nary evolution theories and to investigate the last stages 
of single star evolution as well. We note that only 10 
double white dwarfs were detected during the last 10 
years (Maxted et al. 2000) and about 40 pairs of a white 
dwarf and a neutron star - during 25 years (Thorsett & 
Chakrabarty 1999). We further note that by means of 
gravitational lensing it is possible to detect compact ob- 
ject binaries at the same rate (see Section 5). Moreover, 
this may help to discover binaries containing neutron stars 
with a low magnetic field (not pulsars). Study of these sys- 
tems in combination with binary star evolution theories 
gives a possibility to test the detailed cooling models for 
white dwarfs and neutron stars and the equation of state 
of the latter (Nelemans et al. 2001; Yakovlev et al. 1999). 
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Studying compact object pairs is very important to 
solve a number of astrophysical problems. We present be- 
low several such examples. 

The "Standard candle" of modern cosmology, type 
la supernovae, apparently arises from merging of double 
CO white dwarfs (Webbink 1984; Iben & Tutukov 1984), 
which have not been detected at this time. 

Close double compact objects have to contribute a sig- 
nificant part of the gravitational wave signal at low fre- 
quencies. Thus, white dwarf pairs may be a source of 
the unresolved noise (Evans et al. 1987; Grischuk et al. 
2001). Statistical properties of the closest double compact 
objects in principle determine parameters of the gravita- 
tional wave signal. 

Binary systems consisting of a white dwarf and 
a neutron star (a pulsar) are unique laboratories for 
high-precision tests of general relativity. To date, post- 
Keplerian general relativistic parameters have been mea- 
sured in four such pairs (Thorsett & Chakrabarty 1999). 
This problem may be solved very easily for binaries with 
gravitational self-lensing since they are observed nearly 
edge-on. 

It seems that the detection of black holes forming pairs 
with white dwarfs might bo an extraordinarily important 
result of the search for gravitational lensing in compact 
object binaries. Despite the opinion of most researchers, 
in a certain sense black holes have not been discovered so 
far. Only observational data on the behaviour of matter 
close to the event horizon showing its presence may tes- 
tify that a black hole is identified (Damour 2000). There 
is indirect evidence of the presence of black holes in X-ray 
binaries and galactic cores based on the "mass-size" rela- 
tion expected for them (Cherepashchuk 2001). The hori- 
zon neighbourhood is seen neither in X-ray binaries nor in 
the cores of active galaxies because they are screened out 
by the accreted gas (the accretion rates are very high). 
This means that only black holes accreting usual inter- 
stellar plasma at low rates of lO""*^^ — lO~^^M0/ycar can 
be recognized as objects without a surface and with an 
event horizon (Shvartsman 1971). The black hole compan- 
ions of white dwarfs in the binaries detected by means of 
gravitational lensing coidd be the best objects for horizon 
study and tests of general relativity in the strong field limit 
(Damour 2000). It is very easy to measure the mass and 
size of such a black hole duo to binary cclge-on orientation 
and investigate the radiation of gas near the horizon. This 
matter is accreted from interstellar medium only because 
its transfer from the white dwarf is absent. 

The gravitational lensing effects in binary systems con- 
sisting of compact components were studied previously in 
several papers (Macder 1973; Gould 1995; Qin et al. 1997). 
However, the authors did not go further than estimation 
of the probability of detecting the effect of orientation of 
the binary system (often underestimating it by a factor 
of 2-3), which turned out to be very low. Their general 
conclusion is that the effect cannot be actually recorded. 
Nevertheless, the data accumulated by now on the evolu- 
tion of binary systems and their parameters make it pos- 



sible to define with higher accuracy the probability of de- 
tection of the brightness enhancement in binary systems, 
to find the expected number of such flares and to propose 
in the final analysis the strategy for their search based on 
present-day facilities. Our paper is devoted to performing 
these tasks. 

We examine in Section 2 the distinguishing features of 
light magnification in the systems consisting of two white 
dwarfs, a white dwarf and a neutron star, and also a white 
dwarf and a black hole. In Section 3 we derive probabili- 
ties of recording the effect for its different amplitudes in 
all three cases. In Section 4 the numbers of systems of 
different types which may be detected by means of gravi- 
tational lensing are estimated. In Section 5 we discuss the 
possibilities of the quest for brightness magnification in 
such systems with the aid of the telescope and equipment 
used in the survey SDSS (York et al. 2000). 

2. Light curves in gravitational lensing 

We will consider three types of binary systems consisting 
of 

a) two white dwarfs with masses of 0.7Mq (sometimes 
0.5Mq), 

b) a white dwarf {Mwd = 0.7 Mq) and a neutron star 
{Mns = 1.4M0), 

c) a white dwarf {Mwd = 0.7 Mq) and a black hole 
{Mbh = IOMq). 

It is clear that both components of a pair play the part of 
a gravitational lens alternately. However, we will mention 
at once that the source of radiation will always be the 
white dwarf to which the approximation of a uniformly 
luminous disk can be applied. The amplitude, the shape 
and the duration of flares of the lensed object are deflned 
by the relationships between the parameters of the binary 
system, the major semiaxis a, the masses and radii of the 
components and the orbital plane orientation. 

The luminosity magnification of a uniformly radiating 
disk in gravitational lensing was analysed in detail in the 
papers by Refsdal (1964), Liebes (1964), Byalko (1969), 
Ingel (1972, 1974), Maeder (1973) and Agol (2002). On 
accoimt of the conservation of the surface brightness, the 
task is reduced in the flnal analysis to a study of varia- 
tions of the source image area with allowance made for 
its occultation by an opaque lens, as the source and lens 
are moving with respect to the observer (in our case - 
the orbital motion around the mass centre of the binary 
system). 

Let ri be the distance from the observer to the gravita- 
tional lens, a is the binary system orbit semiaxis (for the 
sake of simplicity we consider only circular motions) , Mi 
is the lens mass, Vg = 2GMi/c^ is its gravitational radius, 
Ri is the lens radius, Rs is the source radius, Mg is the 
source mass. 

In our case for the light path we can use with high 
accuracy an approximation of a broken line, the angle be- 
tween two sections of which is inversely proportional to 
the impact parameter. 
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The gravitational lens equation for a point-like source 
(see, for instance, Zakharov 1997) is 



- ri + CL-p -ai-E\ 
ri 



(1) 



where rj and ^ are the vectors that characterize the lo- 
cation of the points of intersection of the Ught ray with 
the planes of the source and the lens, respectively, 0{S,) is 
the angle of light deflection (see Fig. 1). In our designa- 
tions 6{^) = ^ • |. The main parameter of the task is 
the Einstein- Hvol'son radius r. = «/2ro-^7*-, the radius 

of the ring - the image of the point source, lying on the 
"observer - lens" axis in the lens plane. In our case ri ^ a 
and 



re = \/2rga 

or, in solar units. 



(2) 



Equation (1) has two solutions and hence the source im- 
ages are two points in the lens plane with coordinates: 



a,2 



1±\1 + 



4r| 

d2 



(3) 



where d = \rj\ is the "lens-source" vector projection onto 
the lens plane, i.e. displacement of the source off the 
"observer lens" axis. 

For the magnification coefficient K, i.e. the ratio of 
the sum of brightness of the two images to the point-like 
source brightness, we have (Zakharov 1997) 



K = 



■qdr] 



+ 



r]dri 




To evaluate the possibility of using particular approx- 
imations, it is necessary to estimate the values of the task 
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Fig. 1. Schematic view of gravitational lensing. 



parameters. As has already been mentioned, we consider 
only the white dwarf as a source of radiation. The role 
of the lens is played by a white dwarf (WD), a neu- 
tron star (NS) or a black hole (BH). At typical veloci- 
ties of the orbital motion, 100-300 km/s, the accretion 
luminosity of the BH is 10'' — 10^ as low as the ther- 
mal luminosity of the WD (Ipser & Price 1982). Assuming 
Mwo = 0.7Mq,Mns = 1AMq,Mbh = WMq,Rns = 
15, Rbh = I'g = we have estimated the rest 

of the parameters of binary systems, which are presented 
in Table 1 in solar units. 

Table 1. Einstein radii for different binary systems 



Pair Ml 


Ri 








Rs 




{Rq) 




[Re) 




(Re) 






a 


= 10 a-- 


= 100 




WD-hWD 0.7 


10-2 


1.0 


■ 10-2 3 • 


10-2 


10-2 


WD+NS 1.4 


2 ■ lO-'^ 


1.1 


• 10-2 3.5 


• 10-2 


10-2 


WD-FBH 10 


4 • 10-^ 


3- 


10-2 9.5 


• 10^2 


10-2 



For the WD radius we used the relationship 
(Nauenberg 1972): 



R = 0.01125i?, 



M 



1.454M( 



M 



1.454M, 



.(5) 



One can see that in a system of two WDs, the Einstein 
radius, the sizes of the lens and the source are close. This 
implies that when considering the light magnification one 
has to take into account the source extension and occulta- 
tion of its image by the companion - the lens, at least, at 
a < lOOi?0, which, as we will sec, is our primary interest. 
In binary systems, containing a NS and a BH, eclipsing 
is absent; the source, however, can be regarded as a point 
one only at a > (10— 2O)i?0, and one can apply expression 
(4) to derive the brightness increase. 

It should be noted that when the binary system is 
viewed edge-on, at the moment of conjunction, a configu- 
ration in which the observer, the lens and the source are 
on one line {d = 0) is realized. It was shown (Liebes 1964) 

^1 + as a result of averaging of 



that in this case K ■ 



expression (4) over the surface of the uniformly radiating 
disk. At d 7^ the task gets essentially more complicated 
and has no analytical solution. The accurate relationships 
for the brightness magnification by a point lens, where el- 
liptical Legendre integrals of the first and second kinds 
are used, were derived by Refsdal (1964), Liebes (1964), 
Byalko (1969), Ingel (1972, 1974) and Maeder (1973). 
There is, however, one more case in which a simple ex- 
pression for K exists - the contact between the disk edge 
and the lens centre: 



K = 



1 + 4(^)2 



arcsin 



1 + 4(^)2 



(6) 
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where j3 = 




200 



Time, s 



Fig. 2. Magnification curve for the pair of white dwarfs. 
The crosshatched circle is the gravitational lensing com- 
ponent. 

The light curve of the componcnt-soTircc is the result 
of variation of magnification with time variations of d dur- 
ing the system rotation period. In the approximation of a 
point source, according to (4) 



1 



4r^ 



(7) 



Near the conjunction (P{t) = dQ + {vtY, where do = asinz 
is the minimum value of d, which depends on the orbit 
inclination z, v is the ralative velocity of the system com- 
ponents. At the moment t = Q, the brightness variation 
reaches a maximum and 



= K{0) = i 



1 + 



(8) 



here/3o = ^. 

It is customary to characterize the duration of the flare 
(brightness increase) (see Zakharov & Sazhin 1998) by the 
time T that is taken by the moving source projection onto 
the lens plane to cross the Eintstein ring. It can be easily 
found from the relation d^(|) = t^o + (^)^ = ''e' t^^'*- 



2re 



T 



1 



PI 



(9) 



In a point approximation K{^) = + = 1.34, and 

T defines the flare duration at a level Amo = 0^32 of the 
source light at a quiescent phase. 



Note that the brightness magniflcation of one of the 
components by a factor of K corresponds to magnification 
of the system brightness (which is the subject of detection 
in observations) Kq times, and 



(10) 



where is the lens-source brightness ratio. For a system 
containing two equal white dwarfs = \, Kq = ^{K + 
1), and Amo = 0^17 with respect to the total system 
brightness off the flare. It is self-evident that if the part 
of the lens is played by a neutron star or a black hole, the 
K is practically equal to Kq. 

Present-day telescopes detect rather small light vari- 
ations during short exposures (see Section 5). We are 
guided by these possibilities when determining the min- 
imum amplitude of detectable flares. Say. if the amplitude 
of a flare Amo = 0^2 [K = 1.2), then the trajectory of 
the source projection onto the lens plane does not cross 
the Eintstein ring at all, and such a flare is excluded from 
the analysis. Therefore, it is necessary to introduce ro.i at 
which K = 1.1. Using (4) we derive 



^!£^2 _r.2_ K^-I + KVk^ 
^d' ~^ ~ 2 

hence at K = l.l,/3o.i = 0.60, do.i = 1.67re. Thus 



ro.i = 1.67— Jl- 



1 



(11) 



(12) 



In the general case, when deriving the light curves, 
equation (1) was solved for each point of the source pro- 
file, its image was constructed in the lens plane, and its 
area was found numerically. Since in gravitational lensing, 
the surface brightness remains constant, the light curve is 
determined as variation of the ratio of the image areas and 
the source itself, depending on the position of the binary 
system components with respect to the observer, which 
varies during the period of rotation. The contribution to 
the system radiation of the lens-white dwarf and occul- 
tation by it of a part of the source image was taken into 
account. Fig. 2 shows the brightness variations of the bi- 
nary system consisting of two white dwarfs with a mass 
of O.7M0, a = 14.2i?0 and the orbit inclination i = (i.e. 
do = 0). The similarity of the sizes of the source, the lens 
and the Einstein ring results in formation of a rather com- 
plex light curve. Similar events were analysed by Marsh 
(2001) and Agol (2002). 

Comparison of the light curves in Fig. 3, plotted from 
the use of the results of direct calculations and with the use 
of point approximation (7), shows that they are virtually 
coincident for a > 30Rq. In this case the light curves can 
be constructed by means of a point approximation. The 
exception is the region near the maximum at < 0.25. 

Consider the effect of gravitational lensing in a binary 
system comprising two equal white dwarfs with Mwd = 
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Fig. 4. Simulated light curves of the pair of white dwarfs. 
Fig. 3. Magnification curves of the pair of white dwarfs for 

different semiaxes a. Dashed lines represent the analytical ^^^^i a step of 2.5 • 10"^ from the upper to the lower curve, 
approximation (7). Solid lines for numerical calculation; ^^^^^ jg displaced for convenience. The flare shape varia- 
fine ones for transparent lens, thick ones for opaque lens, ^ions are due to these changes. 



O.SMq and a = WORq. According to Kepler's third law 
(in solar units): 



Table 2. Parameters of the closest pairs able to give rise 
to a flare of amplitude Am 



T = 2.79a^/^(M; + Ms)~^ hours, 



(13) 



its period T = 82.4 days. Flares which amplitude depends 
on i will be observed twice during this time. For do ai > 
Rs, point approximation (8) can be employed to estimate 



the amplitude. Consequently, at i > ^ 
follows from (2) and (8) 



10- 



= 10- 



it 
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WD+NS 




WD+BH 








ro.i 












(Rq) 


(Re) 


(s) 


(Re) 


(Re) 


(s) 


(Re) 


(r2 


4 


14 


44 


0.92 


2.5 


12 


0.13 




15.4 


23 


170 


3.15 


3.2 


41 


0.44 




62.4 


49 


690 


11 


10 


141 


1.52 



1 + 



2r„ 



< 2.27. 



./i + 4^ 

Thus the amplitude of these flares Kq = ^ [K. 
or Amo < 0.54. 

Prom (2) and (9), considering that v = 
obtain under the same condition (point approximation): 

(14) 



-1) < 1.64 



G{Mi + M,) 



4a 
c 



Ml 



9 • 10-6M, ■ 



4.5-10- 



And for the given system r 

662V1 - 2.2 • 107 . sec. 

In Fig. 4 the light curves of this system for different i, 
plotted from the data of numerical computations, are dis- 
played. The light variations arc expressed in stellar magni- 
tudes. The orbit inclination i varies from zero to 2.5-10"'* 



One more step brings us closer to the real condition 
of searching for flares. Having fixed three levels of bright- 
ness increase. 0™2, 0^5 and 1™0, wc determined by two 
methods the minimum sizes of the systems consisting of 
two white dwarfs and also a white dwarf and a neutron 
star whose light rises to these levels (Table 2). Prom the 

relation K = y 1 + 7^ , amin has been found analytically, 
while Uc has been derived by direct calculations. The light 
curves of flares in these edge-on-viewed pairs, which were 
constructed based on numerical computations, are pre- 
sented in Fig. 5a, b. It can be seen from Table 2 and this 
figure that the analytical estimates of the characteristic 
durations of flares tq.i arc close to the results of direct 
calculations. The flares of a white dwarf coupled with a 
black hole constructed in a similar manner are exhibited in 
Fig. 5c (Table 2 contains only amin for this binary system) , 
here the separation of the components a = Si?©, while the 
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amplitudes 0™2,0™5, 1™0 correspond to the orbital plane 
inclinations 0.023, 0.015 and 0.01. It is precisely this value 
of a that was used in the computation since it is close to 
the minimum size of the pair WD-BH, which "survived" 
in the process of emission of gravitation waves during the 
lifetime of the Galaxy. The number of closer pairs drops 
sharply and the probability of their detection is very low 
(see Section 3). 

^ 1 I ' 1 ' ' — I — . ' 1 ' 1 
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Fig. 5. Simulated light curves of the flares as a result 
of gravitational lensing in different pairs: a) WD-WD, 
b) WD-NS, c) WD-BH. Dashed Hues show flares with 
amplitude 1™0, dotted lines are used for flares with am- 
plitude 0™5 and solid lines for flares with amplitude 0™2. 



Fig. 5a, b, c shows that the characteristic durations of 
flares for pairs of white dwarfs are 50-200 s; for binaries 
containing NS and BH, the light variations are shorter, 
30 130 s and 20 40 s, respectively. Singnificant variations 
of the flare shape are noticeable in all the cases. 

Without dwelling on details, note that the set of pa- 
rameters of these binary systems obtained in observations 
- light curves, periods, amplitudes, duration and shape of 
flares, can make it possible to determine physical charac- 
teristics of the components - their masses, sizes, temper- 
atures (e.g. Cherepashchuk & Bogdanov 1995a, b). This 
problem can be resolved completely for two white dwarfs. 
At the same time, the determination of, at least, the mass 
of the lens component and, therefore, its identification as 
a NS or a BH seems to be self-important for the study of 
properties of relativistic objects and their extreme gravi- 
tation fields. 



3. Probability of detection of flares and their 
average characteristics 

The number of systems that can be found as a result of 
brightness increases in gravitation lensing is determined, 
first of all, by the probability F{Ko, t) of recording a flare 
with an amplitude larger than the level Kq, during the 
observing time t of a sample of systems distributed in 
some manner over parameters (orbital separation, masses 
and component brightness). This means that F(Ko,t) is 
the probability of a flare detection from a certain system 
averaged over such a distribution 



F{Ko,t) = J PiKo,t,uj)f{uj)duj, 



(15) 



where uj is the set of parameters, f{uj) is the function of 
their distribution, f2 is the region in which they are speci- 
fied, P{Ko, t, uj) is the probability of recording of the flare 
in a system with the parameters lo. For the sake of simplic- 
ity, with allowance made for the low accuracy of any esti- 
mates associated with statistical characteristics of binary 
systems with compact companions (Masevich & Tutukov 
1988; Lipunov et al. 1996; Nelemans et al. 2001; Fryer et 
al. 1999, where the errors of different models arc defined), 
we restrict ourselves to an analysis of the distribution of 
the systems only in a. Then P{KQ,t, a) is the probability 
of a joint onset of two events, favourable orientation of the 
orbital plane with the semiaxis a, when the level Kq in the 
fiare is exceeded, and the recording of a result if only one 
such a flare occurs during the time t, f{ui)du = f{a)da. 

The first event is consistent with the geometric prob- 
ability pg which can be easily found. It follows from (4), 
(7), (8) and (11) that the values (I'q < do = ^5^, where 

K = Kq + jj-{Ka — 1), correspond to the flares of the 
source with amplitudes exceeding K. The locus of the 
ends of the normals of the orbital planes satisfying this 
condition is a globe belt of width do and radius a/2, its 
area is Sd = nado. The locus of the ends of the normals of 
any orbit is the surface of a sphere of radius | , its area is 
Sa = na\ Hence pg = ^ = ^ = Finally, from (2) 

Pg = 3 • 10-^a-iMf[(3(K)]'\ 

Note that in determining P{K), both analytical expres- 
sion (11) and the numerical relationship can be applied. 
Fig. 6 shows the relationships between geometric proba- 
bilities and the size of a system consisting of two white 
dwarfs which have been derived by these two methods 
(analytically in a point approximation - the dashed line, 
and numerically - the solid one). The contribution of the 
dwarf lens is taken into account here. Minimum sizes of 
systems capable of giving rise to a flare of amplitude K 
correspond to the values Oc in Table 2. The numerical cal- 
culations in systems WD-NS and WD-BH are virtually 
the same as the analytical estimates, but Pg for WD-NS 
reduces to zero &t a = Uc- 

The probability pt of recording a flare during the time 
t is essentially dependent on the relationship of the period 
of recurrence of flares T (given by (13)), the duration of 
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Fig. 6. Geometrical probability versus semiaxis size for 
the pair of wliite dwarfs. Dasiied lines for analytical ap- 
proximation, solid lines for numerical calculation. 



the flare t and the net continuous time of observing t. 
Let us lay down certain conditions which must be satis- 
fied to make the search for flares efficient enough. First of 
all, it is evident that the duration of an elementary con- 
tinuous exposure tg must be close to that of flare r. On 
the other hand, the interval between the exposures must 
be considerably shorter than r. Implying hy t = '^te the 

e 

duration of the total exposure when observing a specified 
sky region, we obtain for the probability pt- 



Pt 



^ at t<T 
1 at t>T, 



which corresponds to 



Pt = 



0.36-^^ — s/2 — fit a > at, 

1 at a < at, 



(16) 



(16') 



whereat = [0.36t(M;+M,)i/2]2/3 is the semiaxis of binary 
with period T = t and we consider at > ac only since 
Pg = as a < ttc- 
Hence 



P{Ko,t,a) =PgPt 

0, at a < ttc, 

3 • 10-3a-^M//^[/?(/r)]-i at ttc < a < at, 

1.08 • 10-3to-^ '^'(y^^-)'''' at a, < a; 

at < a. 



(17) 



Due to the serious anisotropy of pair distributions on 
the sky, the optimum time of observation of a specified sky 
field will be shown to be greater than one night length (at 
least for WD- WD systems). In this case a simple expres- 
sion (16) for Pt no longer holds since observations on indi- 
vidual nights become statisticaly independent. The prob- 
ability of detecting at least one fiare in n nights is given 



by the Bernoulli distribution: 

n 

Pn = i-{i-Pir = J2^-iy+'cip{. 
1 

where pi = Pt\t=9'^ = 3.24(Mi + Ms)1/2q-3/2 jg ^^le de- 
tection probability in 9^ (the observational night average 
length). 

Thus, expression for P{Ko,n,a) takes the following 
form 

P{Kn,n,a) = PgPn 

'0, at a < ttc 

3 • 10-3a-5M//^[/3(ii')]-i, at ac < a < ai 
3 • 10-^ Ml^^[f3{K)]-^ 

n 

x^(- 1)^+1 (3.24)'C4 
1 

. x(M;-hMs)'/2a-i(3m) atai<a, (17') 

where ai = at\t=9h = 2.2{Mi + Msf/^. 

For randomization corresponding to (15) we must con- 
volve P{Kq, t,a) or P{Kq, n, a) with the pair distribution 
function in a. 

Now we turn to the analysis of this distribution. We 
must clearly see the difference in distribution for binary 
systems, progenitors of pairs of compact objects, and the 
finite distribution of these pairs themselves. In the first 
case the results of observations yield a distribution uni- 
form in logarithm of semiaxis (Abt 1983): /(a) oc a~^(l < 
log(a/i?0) < 6). In the second case the situation is not 
so clear because of the uncertainty in quantitative char- 
acteristics determined in the course of different versions 
of population synthesis and because of the qualitative dis- 
tinctions of these versions (Yungelson et al. 1994; Lipunov 
et al. 1996; Saffer et al. 1998; Portegies Zwart & Yungelson 
1998; Fryer et al. 1999; Nclcmans et al. 2001). For instance, 
in different papers the initial mass functions with different 
indices are used; due to the great uncertainty in the loss of 
matter, the supposed masses of BH progenitors vary from 
25Mq to SOMq. The finite numbers of, say, couples WD- 
BH contain an uncertainty of 5 orders of magnitude (see 
Fryer et al. 1999)! Nevertheless, since our analysis is qual- 
itative, one can be content with the main characteristics 
of samples of objects and rough estimates. 

First of all, it will be noted that couples lose energy and 
angular momentum radiating gravitation waves (Landau 
& Livshits 1983; Grishchuk et al. 2001). Since the rate of 
these losses depends dramatically on a (cx a^^), the stars 
in close pairs stick together and leave the sample. Using 
the expression for gravitational luminosity of a binary 
system in a circular orbit ^ = — • ^' 
(Landau & Livshits 1983) and taking into account that 



GMiM, 
2a 



one can 



the total energy of the system E 
readily derive the value of ar of the initial semiaxis of a 
system that will merge during 10^° years, which is the 
life-time of the Galaxy: 

ar = 2RQ[MiM,{Mi+M,)]^/\ 
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Table 3. The closest pairs that survived during the Ufe- 
time of the Galaxy 



Pair 


Mi+Ms 


Or 


Tr 




(Mq) 


(^0) 


(hours) 


WD+WD 


0.7+0.7 


1.82 


5.77 


WD+NS 


0.7+1.4 


2.40 


7.11 


WD+BH 


0.7+10 


5.88 


12.13 



Wc tabulate in Table 3 the values of and corre- 
sponding Tr for the systems under study. 

The results of population synthesis, which are partially 
confirmed by observations, suggest that the density dcstri- 
bution of pairs of compact objects rises with decreasing a 
approximately as -, reaches a maximum at a = ao ~ ctr, 
after which it decreases rapidly to a l^o (Lipimov et 
al. 1996; Hernanz et al. 1997; Saffer et al. 1998; Portegies 
Zwart & Yungelson 1998; Wellstein & Langer 1999). 

It follows from the same papers that another important 
parameter, the maximum size of a binary system amax^ is 
known to be less than 50Rq at small eccentricities. 

With allowance made for the said above, we have used 
the function of distribution in a in a simple form: 



/(a) =j\geg{a) 



aa( 
7lge 



(18) 



at Go < a < an 



From the condition of normalizing / f{a)da = 1 one 

1 

can easily find 7: 



0.437 = 



gp — 1 
2ao 



+ ln^ 



(19) 



Note at once that the share of systems with a < ao ranges 
from 10 % to 30 % at 3 < ao < 10 and a,nax = 50, 30, 
which is close to the results of population synthesis for 
pairs of different type compact objects (Hernanz et al. 
1997; Saffer et al. 1998; Portegiez Zwart & Yungelson 
1998; Nelemans et al. 2001). 

In our further estimations we will use distribution (18), 
varying the parameter ao in the range 3-R0 + lOi?0, at two 
values of the parameter amax ~ 

3Oi?0 and 5Oi?0. 
Having fixed the distribution /(a), now we have to de- 
termine the lower limits of integrating in (15) for different 
types of binary systems. Referring to Table 2, Table 3 and 
Fig. 6 and allowing for (17) and (18), it can be readily seen 
that the lower limits for systems consisting of two white 
dwarfs and a white dwarf in pair with a neutron star co- 
incide with ac and are 14 and 2.5, respectively. However, 
for pairs consisting of a white dwarf and a black hole inte- 
gration is done over the whole range of possible values of 
a, i.e. from 1 to amax- The values of binary periods corre- 
sponding to Uc are 123, 7.6 and 0.85 hours for WD- WD, 
WD-NS and WD-BH systems respectively. This leads to 



the following expressions for flare detection probability ob- 
tained from (16), (16'), (17), (17') and (18). 

For WD-WD pairs 
when f < 9'* < Tc = 123^* since a\t^Qh < Uc and ao < ac, 

F{Ko,t) = X/(a)l-08- 10-3[/3(i^)]-ito-2 

ac 

X [Mi{Mi + M,)]i/2da = 1.08 • lO'^ • 0.437 
x[/3iK)]-H[Mi{Mi + Ms)]^/^ "Ya-^da (20) 
= 5.4 • 10-^[(3{K)]-H[Mi{Mi + M,)]i/2 



2ao ao 



1 _ 1 



when the overall observation time of a specified sky region 
is n nights. 



F{Ko,n) = 7"V(a)3- 10-3M//'[/3(i^)]-i 

ac 

X E (-l)'+^(3.24)'Q(Mi + M,)ia-^da 
= 3 • 10-^0A3'yMl^^[p{K)]-'^ 
E(-l)'+'(3.24)'C;(Mi+M,)i / a-^da 



= 3-10- 



E(-ir+H3.24)»C;(M; + M,)i3^ 



i=l 

■ _ 3i + l 3i+l 

ac amax 



(20') 



For WD-NS pairs 

when t <Tc = 7.^6, since at < ac < ao 



F{Ko,t) = 1.08- 10-3[/3(jr)]-if[MKM; +M,)]V2 

X 7V(a)a-2rfa = 1.08 • lO'^ 

_i (21) 



ao(ao — 1) 

when 7.''6 < i < 9'*, since ac<at < ao 



F{Ko,t) = 3-lQ-^\p{K)]-^Ml'^ 
+0.36(M/ + Msf'H "'Y a-'^f{a)da 
= 3 • 10-3 



7 a-^/^f{a)da 



'1/2 

+0.36<(Mi + M,)i/2 



M, 



1/2 



^[0.36t{Mi+Ms)^/^-al^^]-2[0.71t^/^{Mi+M^y/'^-ay^] 
ao(ao-l) 



Inao- |ln( 0.36t 
X {Ml + M,) 1/2 ) _ ( 0.36t(M, + M,) 1/2 ) + -i- 



(21a) 
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when the overall observation time of a specified sky region 
is n nights, since Oc = 2.5 < atlt^gh. = 2.8 and at < ao 



1 )l^l/2 



F{Ko,n)=3-10-^[/3{K)]-^M; 



J f{a)a-^/^da 



+ ^(_iy+i(3.24)^C;(M, + M,)i / f{a)a-^da 



3 • 10-3 



x(ao(ao - 1)) 



M, 



1/2 



-2[1.48(/lf, + M,f/^ - aj^] + 3.24n(M, + Msf'^ 

n 

X 1^ 2.2(M,+k)^/^ + E (-1)*+'(3.24)*C;(M, + M,)^ 



■ 3(i-l) 



(3.24(M; + M,)V2)-(i-i) _a- 
E(-l)'+^3•24)'C;(M^+M,)i 

(3.24(Mi + M,)i/2)-(-4)-ao 



3(i-l) 



3i-l 



3i+l _ 3i+l 

(^max ) 



For WD-BH pairs 

when t < Tc = 0.''85, since at < ac < a and < ao where 
ttc = 1 (left distribution edge) 

F{K^,t) = im-10-^[(3{K)]-H[Mi{Mi+M,)Y/^ 
X 7°V(a)a"^da = 1.08- 10-3 



+ ln^ 



t[Mi(M, + M«)]V2 

ao(ao-l) O^'^O - 1 + ^) + 5 - ^) 

when 0.''85 < t < 9'*, since 1 < a* < at|t=9h = 4. 



(22) 



F{Ko,t)= 3-10-3[/3(i^)]-iM, 



/ a-y^f{a)da 



+QM{Mi + MsY'H J a-'^f{a)da 

at 



and for at < ao 



FiKo, t) - 3 • 10-3 l^+i^a^ 



1/2 



0.24t(Mi+M3)^''^-1.42t^''^(Mi+M3)^'''^ + | 



oo(ao — 1) 



+OMt{Mi + M,)i/2 



ao(ao — 1) 



In ao 



-| In [0.36t(M, + M,)i/2] _ 2.8t-2/3(M; + M,)-i/3 



+ 2(aQ O^max) 



or for ao < at 



(22a) 



F(i^o, t) = 3 • 10-3 [l3{K) + In 



n -1 



3/2 



- 2a, 



1/2 



-1.41t-3 (Mi + M^)- 6 



,-1/2 



0.18(M; + Ms)l/2t 



3.8t-|(M;+M,)-i-a-2j}; 



(226) 



when the overall observation time of a specified sky region 

is n nights, for ai = 2.2{Mi + MsY/^ = 4.8 < ao, F{Ka,n) 
coincides with (21') where ac is replaced by 1, and for 
ao < 4.8 

F{Ko,n) = 3 - 10-3M//^[/3(i^:)]-i(//(a)a-i/2da 
+ E(-iy+H3.24)^(M;+M,)*Q / a-^f{a)da) 

2—1 ai 

= 3.10-3[/3(if) (^+ln^ 
^'''^ xMV^(^(iar-2ar + |) 



+2ao - 1.35(Mz + M,)-^/^ 



-E(-l)'+H3.24)^(M;+M,)iC^ 

i=l 

(2.2)-^(M, + M,)-^ 



2 

n 3i+l 



(22') 



We use the value of 0™2 for Kq for all systems since it 
will be shown in Section 5 that flares with such amplitudes 
are detectable with the present-day observing facilities. 

The values of F{Ko,t) for Kq = 1.2(Am = tr2) and 
t = 9^ are presented in Table 4 for a range of ao and 
arnax valucs. Tlicsc rough model estimates may change by 
a factor of 1.5-2 depending on the real pattern of the dis- 
tribution of pairs of compact objects along the semiaxes. 

Fig. 7 represents the behaviour of F{t) = F{Ko,t) (for 
Am = ffi"2) as a function oft and n for the values of ao and 
dmax which give approximately the average values among 
those presented in Table 4 for different system types. 

We have also estimated the average durations and pe- 
riods of reccurence of the events considered and their 
dispersions using (20) - (22'). To do this correctly one 
has to convolve t and T not with the /(a) but with the 
P{KQ,t,a)f{a)/F{Ko,t) according to the share each pa- 
rameter region brings to the probability of detecting a 
flare. These values, formally, depend on the values of t 
or n. However, it is clear that their variations cannot in- 
fluence much the real average characteristic. Actually in 
these calculations we use the expression for P{Ko, t, a) and 
F{Ko, t) in the region where they are linear in t, which as 
we will see in Section 5 corresponds to the optimal strat- 
egy for searching for the flares. 

The results are presented in Table 5 for the value of 
amplitude Am > 0™2. The parameters averaged over the 
range of the used ao (from 3Rq to IORq) are in the first 
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Table 4. Probabilities of detection of flares with Am = (K"2 in 9 hours 



Parameters WD + WD WD + NS WD + BH 
(Re) 



do \ Clmax 




50 






30 






50 




30 






50 




30 




10 


1.4 


• 10" 


-5 


1.6 


• 10" 


-5 


1.7 


• 10"* 


2.2 


• 10" 


-4 


9.8 


•10-* 


1.3 


• 10" 


-3 


8 


1.3 


• 10" 


-5 


1.4 


• 10" 


-5 


2.2 


• 10"* 


2.7 


• 10" 


-4 


1.2 


• 10-3 


1.6 


• 10" 


-3 


7 


1.2 


■ 10" 


-5 


1.3 


■ 10" 


-5 


2.4 


■ 10-* 


3.1 


■ 10^ 


-4 


1.4 


■ 10-3 


1.8 


• 10" 


-3 


5.88 


1.1 


• 10" 


-5 


1.2 


• 10" 


-5 


3.0 


• 10--* 


3.7 


• 10" 


-4 


1.6 


• 10-3 


2.0 


• 10" 


-3 


4 


9.9 


• 10" 


-6 


1.0 


• 10" 


-5 


4.2 


• 10"* 


5.0 


• 10" 


-4 


2.2 


• 10-3 


2.7 


• 10" 


-3 


3 


9.1 


• 10" 


-6 


9.2 


• 10" 


-6 


4.9 


•10-* 


5.9 


• 10" 


-4 


2.6 


• 10-3 


3.1 


• 10" 


-3 


Average 
estimate 






10" 


-5 








4- 


10"* 








2 • 


10-3 







WD-BH 

WD-NS 

WD-WDJ,^ 














')^' 
/' \ 
/' ' 








/ j 

y 


/ 

/ 

/ 






/ 1 








/ ! 

/ , 









10^ 



1 hours 9/1 nights 10 50 



Total field exposure time 

Fig. 7. Probability ¥{€) of detecting a flare with Am > 
CP 2 as a function of the overall observation time of a spec- 
ified sky region measured either in hours or in nights. 

and third columns. Variations of flare duration with vary- 
ing ao do not exceed a factor of 2. In columns 2 and 4 the 
limits of the intervals which contain parameters of 90 % 
of flares are given. The minimum values of the parame- 
ters correspond to the closest systems of all types which 
are capable of giving rise to a flare with the amplitude 
Am > 0™2. In Table 5 we give parameters of events that 
one can hope to discover in observations. 

4. The number of pairs being investigated in the 
Galaxy 

As we have already noted, the present-day numbers of 
different pairs of compact objects are determined in the 

framework of different versions of population synthesis 
(Lipunov et al. 1996; Wellstein & Langer 1999; Fryer et 



al. 1999; Nelemans et al. 2001). However, the main ob- 
jective of this research is the study of systems which arc 
of interest within the scope of gravitational-wave astron- 
omy and of the problem of origin of gamma-ray bursts. 
These are pairs of neutron stars, neutron stars with black 
holes and essentially rarely white dwarfs with black holes. 
By now extensive work has been done on the modeling of 
the origin and evolution of pairs of white dwarfs. It goes 
without saying that these binary systems are best under- 
stood; comparison with observations can be made for them 
(Hernanz et al. 1997; Nelemans et al. 2001). There are at 
least two relatively reliable estimates of the total number 
of binary white dwarfs in the Galaxy, 5 • 10® (Lipunov et 
al. 1996) and 2.5 • 10* (Nelemans et al. 2001). In princi- 
ple, these data are sufficient to determine the number of 
gravitational-lens flares that can be detected in this sam- 
ple. Nevertheless, we have estimated the number of binary 
white dwarfs and also white dwarfs in pairs with neutron 
stars and black holes in the Galaxy in the framework of as- 
sumptions of the progenitors of binary systems consisting 
of compact objects. We have followed Bcthc and Brown 
(1998, 1999) who have determined the number of massive 
binaries containing neutron stars and black holes. 



Table 5. Parameters of flares of amplitude Am > 0™2 



Duration (s) 


Period of 


recurrence (h) 


average 


interval of 


average 


interval of 




90% 




90% 


WD+BH 50 


10-135 


15 


1-48 


WD+NS 60 


20-120 


42 


8-120 


WD-FWD 135 


90-180 


400 


120-500 



We will base our discussion on the following assump- 
tions (as previously, solar units arc used). 

1. The mass distribution of progenitor stars (initial 
mass function) obeys Salpeter's (1955) law in the interval 
O.K M < 120, i.e. 

dW{M) = w{M)dM = 1.35(^)-2-3^dM , (23) 



G.M. Beskin and A.V. Tuntsov: Gravitational lensing in compact binaries 



11 



where dW{M) is the probabihty that the mass of a star 
will be in the interval M ^ M + dM, w{M) is the proba- 
bility density. 

At the observed star formation rate ^j^^ = 
0.9(^)-2-35year-i and the Galaxy age of 10^° years, the 
number of stars is A'' ~ 1.5 • 10^^. 

2. About 100 % of stars are members of binary systems. 
This evaluation is used in the population synthesis as well 
as 50 % (see for instance, Nelemans et al. 2001). 

3. The distribution density in the mass ratio, q = 
Mg/Mi < 1, is constant, i.e. ^ cx const ( e.g. Portegies 
Zwart & Yungelson 1998). 

4. The initial masses of progenitor stars must lie within 
the following intervals to give: 

white dwarf - 1 < M < 10, neutron star - 10 < M < 25, 
black hole - 25 < M. 

In the last case the uncertainty is rather great; the val- 
ues 40Mq (Van den Heuvel & Habets 1984), and 60Mq 
(Woosly et al. 1995) are also used. As for us, we follow 
Portegies Zwart et al. (1997). 

5. The neutron star formed in a supernova explosion 
in a binary system receives a kick. As Cordcs & Chcrnoff 
(1997) have shown, its distribution is well approximated 
by the sum of two Gaussians with standard 175 km/s and 
700 km/s, the share of the former being 80 % while that of 
the latter being 20 %. With these parameters 43 % of bi- 
nary systems survive after the explosion (Bethe & Brown 
1998). As a black hole is formed, one can believe that the 
kick will be lower, at least, inversely proportional to its 
mass (Fryer et al. 1999). For this case, using the relation- 
ships between the probability of decay and the additional 
velocity derived by Bethe & Brown (1998), we have found 
that 87 % of binaries survive after the formation of a black 
hole. 

Now we turn to estimation of the number of pairs of 
different types. 

For the probability dcf) to obtain a binary system with 
the mass of the more massive companion (in our case the 
object is a lens) in the interval Mi ^ Mi + dMi and the 
mass ratio in the interval q~ q + dq, considering (23) and 
the homogeneity of the distribution in q {Mi is measured 
in O.IMq), we have 



# = dqdW = l.ZbiMiy^-^^dqdMi. 
The number of pairs of type i is Ni 



(24) 

hliNhin, i = d 
for a pair of white dwarfs, i = n and i = b for binaries 
containing a white dwarf with a neutron star and a black 
hole, respectively; 7, is the pair survival probability after 
the supernova explosion. 
Clearly 

J Jd<j) = 1.35j J{Mi)-^-^^dqdMi, (25) 



, Mi 



where Qi is the interval where the values of q lie, and M 
is the interval of Mi values. Determine Qi for each type 

of pairs, taking into account that 10 < Mg = qMi < 100 
(assumption 5). Then < g < 1 for a pair WD-WD, 



^ < q < ^ ior & pair WD NS and WD-BH. Thus, 
proceeding from (25), we obtain for a pair of white dwarfs 

-2.35, 

'nv±i J 

10 

(26) 



<j,d = 1.35 / {Mi)-^-^^dMi J dq 
10 20, 



= 1.35/(M0-2-3^(l-^) dMi. 
10 V '/ 

Hence 



100 100 
= 1.35 / Mj-'^-^^dMi - 13.35 / M^^-^^dMi 

10 10 
= 0.043 - 0.025 = 0.018, 

and allowing for the fact that the number of pairs in the 
Galaxy is N^in = 0.5N, obtain for the number of binary 
white dwarfe Nd = 1.35 • 10^ pairs. 

For white dwarfs in pair with neutron stars 



250 



1.35 / M-'^-^^dMi J dq 



100 

250 



100 

and finally 

90-1.35 



250 

= 1.35 / M-'^-^^j^dMi = 90 • 1.35 / M'^-^^dMi, 
' 100 



2.35 





250" 








100_ 



51.7-1.76-10-^ = 9.1-10"^. 



Because 57 % of binary systems decay after the supernova 
explosion which results in formation of a neutron star, 
N„ ~ 9.1 • lO"'^ • 0.43 • 0.5A^ ~ 2.9 • 10^ pairs. 
The same is for pairs containing black holes: 



= 51.7 





1200" 


M-2.35 






250 



= 51.7- 2.26- 10"^ = 1.2- 10" 



Since the probability of their survival is 87%, then A^^ 
1.2 • 10-^ • 0.87 • 0.5 • A/" ~ 7.8 • 10^ pairs. 

It should be emphasized that the above-given esti- 
mates are defined in many respects by the choice of the 
power index k in the initial function of stellar masses. 
Clearly they decrease with increasing k. It is likely that 
the version with k = 2.35 of Salpeter (1955) is most con- 
sistent with current observational data (see, for instance, 
Grishchuk et al. 2001; Raguzova 2001), but, nevertheless, 
fc = 2.5 and also k = 2.7 are used in the population syn- 
thesis as well (Bethe & Brown 1998, 1999; Nelemans et al. 
2001). 

To make the picture complete, we present in Table 6 
the estimates of the numbers of pairs in the Galaxy for 
the mentioned k. 

On the other hand, in order to determine the accu- 
racy of our estimates, one can compare them (for binary 
white dwarfs) with the results of population synthesis. In 
particular, the number of pairs of this type is 4.5 • 10* 
(Lipunov et al. 1996) and 2.5-10* (Nelemans et al. 2001). 
Note that the last number was obtained for k = 2.7 and it 
is 2 times as low as our rough value from Table 6. The 2- 
4-time differences are reasonable, taking into account the 
uncertainty in the results of modeling of evolution and the 
qualitative character of our estimates. 
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Table 6. Estimates of the number of pairs of compact 
objects for different mass functions 



k 


WD- 


-WD 


WD-NS 


WD-BH 


2.35 


1.3 ■ 


10^ 


2.7- lO'^ 


8-10® 


2.5 


8.7- 


10* 


1.6 • 10^ 


4 • 10"^ 


2.7 


5.5- 


10* 


6-10® 


1.2 • 10® 



Unfortunately, there is no data to compare with for 
binaries other than WD- WD. We suggest only that our 
estimation range may be 2 — 3 times narrower. Therefore, 
the real number of objects may be a few times higher. 

Table 7. The local densities of objects derived from Table 
6 data and distirbution (27) (in pc~^). 



k WD-WD WD-NS WD-BH 



5. Strategy of search and the number of 
detectable objects. Conclusions 

To determine the number of objects that may be detected 
using a certain telescope one has to know not only the 
overall number but also their distribution in the Galaxy, 
their luminosity function and the distribution of absorbing 
matter as well. 

We used the star distribution in the Galaxy of the 
following form: 



-fi. _i£i 



(27) 



where R and z are the point cylindrical coordinates in 
the Galaxy, H and h are radial and vertical distance 
scales, respectively, for which we used the following val- 
ues: H = 8 kpc and h = 250 pc (sec Dehnen & Binney 
1998; Nelemans et al. 2001). Adopting the Sun's coordi- 
nates i?0 = 8.5 kpc and Zq = 30 pc we have: 

Po = Pioce~^e~fr. 

The overall number of objects in the Galaxy is ex- 
pressed by the integral 

oo oo 

N= dz 2TrRpoe~^e~^dR, 



-oo 



which yields 



Plo 



and 



We H 



N 



6.56-1011 



pc 



p{R,z) = 



N 



1^1- 



6.56- 10" 



pc 



(28) 



(29) 



We used distributions for WD-NS and WD-BH in the 
same form with the only difference in normalization com- 
ing from their fewer overall number in the Milky Way. The 
local densities corresponding to the overall numbers given 
in Table 6 are presented in Table 7. 

Note that the part of WD-WD pairs among the com- 
plete number of WDs is (7-25) % (Iben et al. 1997; Fryer 
et al. 1999) and the WD density in the Solar neighbour- 
hood is (4 - 20) • 10~^ pc~^ (sec Nclcmans et al. 2001 and 
references therein). This gives us a range of local density 
of 3 • 10~^ — 5 • 10~^ pc~^. It is very close to our estimate 
shown in Table 7 and the most optimistic values of density 
may be 2 — 2.5 times more than those. 



2.35 2.0 • 10"^ 
2.5 1.3 • 10"^ 
2.7 8.4 • 10"" 



4.1 • IQ-^ 1.2 • 10"^ 
2.4 • IQ-^ 6.1 • 10"® 
9.1 • 10"® 1.8 • 10"® 



The insterstcllar absorption is proportional to the in- 
terstellar medium density which is distributed in the 
Galaxy by the same (27) law with ft, = 80 pc (Dehnen 
& Binney 1998). The normalization was done according 
to the local average value of absorption Ag, = l^lGkpc^^ 
in the g' band of the SDSS equipment (Schlegel et al. 
1998). 

The WD luminosity function is well known and is ex- 
plained fairly well in the frame of their cooling theory. We 
used the luminosity function of Oswalt et al. (1996). It 
is readily seen from this that more than half of the WDs 
are brighter than 13™5 and only less than 20 percent are 
fainter than 15™, where the function reaches its maximum. 
Thus, the magnitudes of the majority of WDs do not suf- 
fer much from the bolometrical corrections which is a few 
tenth for the stars of these spectral classes (see, e.g. Allen 
1973). 

In other binary systems, the situation is more compli- 
cated. It is only clear that the optical emission from both 
a NS and a BH can be associated with interstellar gas 
accretion. In the former case it is likely to be negligibly 
small because of the high orbital velocity of the NS (> 100 
km/s) (Shapiro & Teukolsky 1983). At the same time the 
orbital velocity of a BH in a pair with a WD is about 50 
km/s, and the luminosity of accretion plasma can reach 
Iq30 gj-g/g by moderate-optimistic estimates (Shvartsman 
1971; Beskin & Karpov 2002) and will be by an order of 
magnitude less by a pessimistic one (Ipser & Price 1982). 
Apparently, its contribution to the total luminosity of the 
system is insufficient to correct the depth of the space of 
detection. However, this level of optical emission makes 
possible searching for its fast variations with a microsec- 
ond time resolution for investigation of accretion processes 
and strong gravitational fields near the horizon of events 
(Shvartsman 1971; Beskin et al. 1997; Beskin & Karpov 
2002). 

Now we have everything we need to calculate the on- 
sky density of objects in a given direction {l,b): 



ail,b) 



cos bdldb 

= Jr'^p{R,z)I{M)dr 



(30) 



where I{M) is a cumulative luminosity function from 
Oswald et al. (1996), M = - 51g(r/10) - A{r,l,b), 
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miim is an equipment limiting apparent stellar magnitude, 
A{r,l,b) is a total extinction in the {l,b) direction up to 
the distance r, and 

^; = 00 + r sin b, 



0.8 



R = ^ Rq + cos^ I — 2i?0r cos b cos I 

are the point cylindrical coordinates expressed in terms of 
its galactic coordinates (Z, b) and distance r. 

Unfortunately, this integration cannot be done analyt- 
ically and so we will turn at once to a particular version of 
the equipment being discussed. In our view, an ideal tool 
for carrying out the proposed programme is currently the 
2.5-metre telescope at APO (New Mexico). It has been 
used to accomplish one of the most promising projects - 
the Sloan Digital Sky Survey (SDSS). The telescope has 
a field of 3°, in which 30 CCD chips of 2048 x 2048 pixels 
are mounted (York et al. 2000). Sky scanning at a speed of 
its movement along a strip 2? 5 wide is performed within 
the framework of the project. During an exposure of 55 s, 
the limiting stellar magnitude in the g' filter (it is close to 
the B filter of the Johnson system) is 23™2 at a signal-to- 
noise ratio of ~ 5 (Oke & Gunn 1983; York et al. 2000). 
We have constructed a relationship between amplitude of 
flares recorded at a 5a level and stellar magnitude of a 
quiet object in the g' band (Fig. 8). In so doing, we used 
the parameters corresponding to the SDSS project: the 
seeing - 0.8 arcsec, the quantum efficiency of the CCD 
matrix - 80 %. 

By comparing the data of Fig. 8 with the mean param- 
eters of flares (Table 5) and their light curves for different 
binary systems (Fig. 5), it is clear that with a particular 
exposure of time ^ 10 — 30 s, the flares with Am > 0™2 
of any pairs that we have examined can be recorded at a 
significance level of 5 • 10~^ — 10~^ (duration of any flare 
is a few te)- In fact, the level of the flare detection will be 
better due to registration in the SDSS of the fleld in sev- 
eral colour bands [u' ,g' ,r' ,z') and the limiting r' stellar 
magnitude is also 23™ Moreover, the use and comparison 
of the data in the g' and r' fllters give the possibility of 
easy separation of cosmic ray traces. 

On the basis of (30) we have built a map of the ex- 
pected on-sky densities of WD- WD pairs up to mum = 
23™ which is presented in Fig. 9. The value of the local 
density corresponding to the Salpeter mass-spectrum in- 
dex of 2.35 is used. One can see from the map that the 
maximal density regions have galactic latitudes of 5 — 10 
degrees and a question might arise if it is possible to ob- 
serve in such star-rich regions of the sky. However, even 
in these latitudes the average distance between stars of up 
to 23™ is about 3-4 arcseconds (sec e.g., Zombeck, 1982, 
p. 34) and thus, they are resolved with the SDSS seeing 
and this crowding therefore cannot have much negative 
impact on observations with this equipment. 

To derive the expected number of events of the dis- 
cussed type Nexp one multiplies the probability of detect- 
ing it from a randomly chosen pair F{t) for Am = 0^"2 
(see (15), (20)-(22') and Fig.7) by the number D{n) of 



. T-H 



i-i 



Q 



0.6 



0.4 



0.2 



Exposure time 10 seconds 
Exposure 1 00 seconds 
Exposure time 30 seconds 



18 20 22 24 

Object stellar magnitude g' 

Fig. 8. Magnitudes of detectable flares in the SDSS versus 
object brightness 



.50, . 








-IQ:. 



500 



Fig. 9. Simulated on-sky distribution of WD- WD pairs, 
fniim = 23™. Maximum density value is 490 pairs per 
square degree. 

such pairs observed during a certain observational pro- 
gramme that covered the celestial sphere part O: 



exp - F{t) X £)(0), 



where 



D{Q) = j a{l,b)cosbdldb . 



(31) 



(32) 



It is clear, however, that when applying to reality one 
deals with a telescope of limited field of view and the ob- 
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servational time provided for a programme is limited as 
well, so increasing the total exposure of a field to increase 
F{t) leads to decreasing the number of such fields in the 
programme. It is also evident that one has to monitor at 
first the best, i.e. maximal density a{l, b) fields first. Fig. 
10 represents the number of pairs of each type D{Q), nu- 
merically calculated from (32), in fields with a size of 
2? 5 X 2? 5 (field of view of SDSS) after sorting them in 
decending order of pair density. 
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Fig. 10. Expected number of objects D(ri) in the best 
(with the highest density of pairs) 2? 5 x 2? 5 fields. 

Now let the programme last for 5 years. Considering 
that the matter in question is detection of rather faint ob- 
jects, observations have to be made on dark nights only 
at moon phases within iS'^ from the new moon, which 
gives approximately 36% of night time. Assuming the av- 
erage observational night of 9'' the 5-year resource of time 
then will make Uohs = 657 nights or Tohs = 5913 hours. 
Thus, as every field is observed for t hours or n nights (and 
J7 = Tobs/t, n-obs/n) the expected number of detected pairs 



or 



D 



Tohs 

nobs 



F(n). 



(33) 



(33') 



These quantities derived from Figs. 7 and 10 are pre- 
sented in Fig. 11. One can see from this figure that the 
optimal total exposure per field for WD-WD pairs is 6 — 7 
nights, and it is about 1 night when searching for WD-NS 
or WD-BH systems. 
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Fig. 11. Expected detection numbers in the 5-year pro- 
gramme as a function of total exposure per field. 



The expected numbers of detections with optimal 
strategies along with the optimal total exposures per field 
are given at Table 8. 



Table 8. The number of pairs of compact objects de- 
tectable over 5 years and the optimal total exposures per 
field 



k 


WD-WD 


WD-NS 


WD-BH 


2.35 


22 


9 


16 


2.5 


15 


5 


8 


2.7 


9 


2 


3 


Optimal 
exposure time 
per 2? 5 X 2? 5 field 


6-7 nights 


1-2 nights 


6-9 hours 



A new project entitled "The Dark Matter Telescope" 
(tittp://www. dmtelescope.org/index. htm) being de- 
veloped involved 8.4 (effectively 6.9) metre telescope 
with a field of view slightly wider than that of the 
SDSS. The use of it in a manner similar to that de- 
scribed above increases the number of detectable ob- 
jects by a factor of 2.2. Other wide field telescopes 
of 4 metres in diameter being built now, LAMOST 



( littp : //www, greenwich-observatory . co .uk/lamost . 
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html) and VISTA (tLttp://www. vista. ac.uk), could de- 
tect 1.5 times more binaries than the SDSS. 

Thus, for searching for the flares caused by gravitation 
lensing in binary systems with compact companions, the 
facilities of the SDSS can be used, having reduced the time 
of an individual exposure to 10-20 s. It is obvious that after 
recording a flare in a certain object, another telescope of 
similar class must be involved in observations. This object 
should be monitored to prove the effect and for detailed 
investigation of the detected system. 

Thus, within the frame of the proposed programme, 
it will be possible to find several dozen binary systems 
comprised of compact objects. It will be recalled that in 
numerous surveys of microlensing, the number of reliably 
detected events amounts to about 350 (Alcock 1999), only 
in two cases the matter in question being of massive lenses 
(black holes?) (Bennett et al. 2001). When applying the 
technique suggested, there is a possibility of detecting with 
absolute assurance more than 10 black holes with open 
event horizons. An investigation of these objects with high 
time resolution within the MANIA experiment (Beskin et 
al. 1997) may permit at last the detection of observational 
evidence of extreme gravitational fields. 
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